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ABSTRACT 

We describe a general yet simple method to analyse the propagation of nuclear reaction rate 
uncertainties in a stellar nucleosynthesis and mixing context. The method combines post- 
processing nucleosynthesis and mixing calculations with a Monte Carlo scheme. With this 
approach we reanalyze the dependence of theoretical oxygen isotopic ratio predictions in first 
dredge-up red giant branch stars in a systematic way. Such predictions are important to the 
interpretation of pre-solar AI2O3 grains fro m meteorites. The rea ction rates with uncertainties 
were taken from the NACRE compilation (4ngulo et al.| 1999). We include seven reaction 



rates in our systematic analysis of stellar models with initial masses from 1 to 3 Mq. We find 
that the uncertainty of the reaction rate for reaction ^*0(p, q;)^^N typically causes an error in 
the theoretical ^'^O/^^O ratio of ~ +20/-5 per cent. The error of the ^'^O/^'^'O prediction 
is 10-40 per cent depending on the stellar mass, and is persistently dominated by the com- 
paratively small uncertainty of the ^^0(p, 7)^^F reaction. With the new estimates on reac- 
tion rate uncertainties by the NACRE compilation, the p-capture reactions ^^0(p, a)^"'N and 
^^0(p, 7)^*F have virtually no impact on theoretical predictions for stellar mass ^ 1.5 Mq. 
However, the uncertainty in ^^0(p, a)^^N has an effect comparable to or greater than that 
of ^^0(p, 7)^^F for masses > 1.5 Mq, where core mixing and subsequent envelope mixing 
interact. In these cases where core mixing complicates post-dredge-up surface abundances, 
uncertainty in other reactions have a secondary but noticeable effect on surface abundances. 

Key words: abundances - stars, nucleosynthesis 



. ;-( 1 INTRODUCTION 



X 



Observations of abundances and abundance ratios in stars yield 
powerful constraints on models of internal processes in stars. In 
particular, the predicted abundance evolution of stars of a given 
mass and metallicity depends on both mixing and nuclear burning 
processes. While it is customary to provide observational results 
with some estimate or analysis of the associated errors, this is sel- 
domly done for theoretical abundance predictions. This is partly 
due to the fact that quantitative uncertainties on the input physics is 
not available and that theoretical model uncertainties (in particular 
for poorly understood mixing processes) are hard to make. Finally, 
a systematic error propagation in stellar models can be computa- 
tionally demanding and the effort has often not been justified in 
the past. However, new spectroscopic observations with unprece- 
dented precision will be available in the near future, and one has 
the impression that stellar physics in general is entering a new 
high-precision era. In this context it seems necessary to reconsider 
the question of abundance prediction uncertainties in a quantitative 
way. 
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As an initial step in this direction, we have chosen specifically 
the dependence of oxygen isotopic ratio predictions in giant stars 
on nuclear reaction rate uncertainties. The oxygen isotopic ratios 
are of particular interest because they are affected not only by nu- 
cleosynthesis but also by mixing processes which are not yet very 
well understood. The latter is evident from the spectroscopic prop- 
erties of giant stars whic h are not entirely consistent with standard 
stellar evolution models ([Harris, Lambert. & Smitl] |l985|; Harris et 



1987| ; |Boothroyd & Sackmann] |199S| ). Therefore it is important 



al. 

to quantitatively know the uncertainties arising from nucleosynthe- 
sis in order to help constrain mixing processes using observations. 
In addition, the oxygen isotope predictions are important for the 
interpretation of pre-solar corundum (AI2O3) from meteorites for 



Huss et. al 


1994; 


Nitfler et al. 


1997) 



(e.g. 

During the main sequence (MS) evolution of low- and 
intermediate-mass stars, partial H-burning in the envelope produces 
"O, ^^C, and "N while ^^C, ^^N, and ^*0 are destroyed (Fig.[l]). 
The abundance of ^^O in the envelope is essentially unchanged for 
stellar mass ^ 1.5 Mq. For more massive stars, core convection 
on the early MS reaches out to ~ 25 per cent of the stellar mass, 
leaving behind core processed material as it retreats, including de- 
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stroyed ^''O (e.g. [Schaller, Schaerer, Meynet, & Maeder 1992 ). Ev- 
idence of the retreating core convection is seen in tiie 2 and 3 Mq 
cases in Fig. |l| at tiie bottom of tiie envelope where the H and ^''O 
are depleted. 

At the end of the MS phase hydrogen is exhausted in the core, 
causing the core to contract and the envelope to expand and cool. 
Envelope convection descends into the star as the envelope temper- 
ature drops. Regions previously affected by cool H-burning are ho- 
mogenised, leading to a change of surface abundances (1°* dredge- 
up, Idup hereafter). As a result of this mixing, theory predicts the 
i6Qyi7Q decrease from initially ~ 2600 to a few hundred, 
depending on stellar mass while the ^^^0/^*0 ratio increases only 
marginally by ^ 20 per cent from the initial value of about 500, 
weakly dependent on mass. For many stars, this is in fair agreement 



with s pectroscopic observations. For example, Harris & Lambert 
( |1984 find for K-giants ^'^O/^^O ratios in the range 300 .. . 1000 
andT^O/^^O ~ 425 . . . 600. Similarly, some corundum oxygen 
isotopes are in rough agre ement with standard Idup r nodel pre- 



dictions (e.g. SEAL203 in Choi, Huss, & Wasserburg (1998) has 



0/"0 = 355 and a solar ^'^O/^^O ratio) while many others 



require either non-solar initial 'O/ O rat ios (Huss et. al 



Weaver 



1994 



Boothrov d. Sackmann. & Wasserburg 1994|^Timmes, Woosley, & 



1995| ) or non-standard mixing processes 



The dependence of oxygen isotopic ratios on t he ^^O proton 
capture nuclear reactions has been studied before by El Eid| ( 1994) 



and Boothroyd, Sackmann, & Wasserburg (1994). The y considered 
t he eff ect of a new determination of these reactions by L 



^ andre et al 

( 1989| ) and found that the uncertainty of the ^^O/^^O predictions 



is dominated by the uncertainty of these rates for stars with initial 
masses larger than ~ 2.5Mq. 

In this paper we include seven reactions in a systematic anal- 
ysis of the influence of rate uncertainties on the predicted oxygen 
isotopic ratios in red giant branch (RGB) stars. In addition to ad- 
dressing this scientific question, we also want to demonstrate and 
test the method of evaluating model error propagation described 
here. Sect. 2 describes the method and computations. In Sect. 3 the 
uncertainty in the modeled oxygen isotopic ratios of RGB stars is 
analysed for significant reaction rate uncertainties and comparisons 
are made to observations and other models. A summary is in the fi- 
nal section. 



network, including the PP-chains and the CNO tri-cycle, is incor- 
porated. 

The treatment of mixing, in particular the use of convective 
overshooting, is not essential to this study because we are inter- 
ested here in the differential impact of the uncertainties of nuclear 
reaction rates. However, there is some agreement that the hydrody- 
namical properties of convection inevitably result in some turbulent 
mixing into the neighboring stable layers, defined by the convective 
boundary where the acceleration of streams disappears. Models of 
convection in the hydrodynamical framework predict that the turbu- 



lent velocity field decays roughly exponentially (e.g. Xiong 1985 
Freytag, Ludwig, & Steffen|l99( ; Asida & Arnett 2000 ). However, 
the extent of overshooting is not the same at all convective bound- 
aries. This is evident from hydrodynamica l simulations, for exam- 
ple by iFreytag, Ludwig, & Stefferl (1 19961) who find fov = 0.25 



for the shallo w surface convection of A-stars and fpy = 1.0 for 
white dwarfs, 
found that fov 



Herwig, Bloecker, Schoenbemer, & El Fid (1997) 



0. 02 reproduces the results of Schaller, Schaerer, 



Meynet, & Maeder (1992) who fitted the observed width of the 
main sequence with an instantaneous treatment of overshooting of 
0.2Hp. This order of magnitude for core convect ion overshootin g 
was confirmed by 2D hydrodynamical models by Deupree ( 2000 ). 
For a situation with an extended, deep envelope, like during the first 
dredge-up, an estimate of overshoot i ng effi ciency is more difficult. 
Alongi, Bertelli, Bressan, & Chiosi (1991) found that the position 
of the RGB luminosity bump can be aligned with observation with 
an envelope overshooting th ree times the value requi r ed for the core 



oversh ooting (see, however, Cassisi & Salaris 



1997). Blocker etal 



(1998) found that overshooting at the bottom of the solar convec- 
tion zone in excess o f fpy = 0.07 would ca use too much lithium 

1999). In view of the am- 



Schlattl & Weiss 



destruction (see also 
biguous evidence we decided to use the same value as for the main 
sequence core convection at the bottom of the descending envelope 
convection during the first dredge-up evolution. 

For the core convection we use no overshooting in the 1 Mq 
cases, fov ~ 0.008 in the 1.5 Mq sequence and fov = 0.016 for 
the 2 Mq and 3 Mq cases. The mixing efficiency in the envelope 
is chosen t o be fenv = 0. 016 for all cases. The Reimer's mass loss 
formulae (Reimers 1975) with -q — 0.5 for M < 1.7 Mq and 
rj — 1.0 for M > 1.7 Mq is employed, starting at the bottom of 
the RGB. 



2 METHOD AND MODELS 

The nucleosynthesis and mixing considered here do not alter the 
structural evolution. Therefore we can efficiently post-process a 
time sequence of full stellar models to analyse the behaviour of 
nucleosynthesis and mixing in the envelope under variations of the 
input nuclear reaction rates. A sequence is post-processed many 
times, each time with a different set of values for the reaction rates 
(within published uncertainty) to obtain a Monte Carlo (MC) esti- 
mate of error propagation into the surface abundance predictions. 
In the following subsections we describe the stages of this method 
in more detail. 



2.1 Evolution code 

A ti me sequence of stellar model s is generated by the EVOL code 
(see Herwig 200C; Blocker 1995, and references therein). At each 



time step the four partial differential equations of stellar structure 
are solved on a ID Lagrang i an gri d, with the OPAL opacity ta- 
bles from ^glesias & Rogers) ( 1996| ). A standard nuclear reaction 



2.2 Post-processing code and the Monte Carlo Sclieme 

Our post-processing code takes the stellar structure from the full 
stellar evolution model sequences generated by EVOL as input and 
recomputes abundance changes due to nucleosynthesis and mixing. 
The code is fully implicit, iterative, and couples mixing with nucle- 
osynthesis using adaptive step size and full error control. 

All reactions in the po st-processing code use NACRE reaction 
rates (Angulo et al, 1999), which include estimates on upper and 
lower limits. We have used the electronic web toolWe/gen^^Jorissen 



& Goriely |2001| ) to retrieve the tabulated rates. The full network of 
nuclear reactions is the same as in the full stellar evolution code, 
but rates for a predefined set of seven reactions were systematically 
altered. 

If one is interested in knowing the abundance uncertainty of a 
particular species, then inspection of its production and destruction 
reactions may reveal to which reaction uncertainty the abundance 
is most sensitive. For example, ^**0 is mainly destroyed rather than 
produced in the envelope. Hence, the surface abundance of af- 
ter Idup ought to be most sensitive to the one of the two major 



© 2002 RAS, MNRAS 000, 



Oxygen isotopic ratios in first dredge-up red giant stars and nuclear reaction rate uncertainties revisited 



reactions wliich consume it. Reaction ^*0(p,a)^^N has a larger 
nominal rate than ^^0(p, 7)^^F and accordingly this reaction must 
be most important to ^*0. While such a qualitative consideration 
is straightforward it is harder to estimate the quantitative effect 
of the uncertainty of this reaction because the combined effects 
of ^**0(p, a) and the other reactions when the temperature profile 
evolves before equilibrium is reached and core convection on the 
early main sequence may affect the envelope abundances. 

Manipulating the rates of one reaction that is part of a larger 
network will not necessarily reveal the whole picture, especially 
when equilibrium abundances are not achieved. This is true in par- 
ticular for a general network situation in which reaction flows may 
be redirected due to uncertainties of reactions that are involved in 
branchings. In order to have a tool available which can readily be 
applied to any nucleosynthesis and mixing scenario we decided to 
take a general and flexible approach by combining post-process cal- 
culations with an MC scheme, as described below. With this ap- 
proach we could, for example, also study more complex situations 
like the dependence of the s-process in asymptotic giant branch 
(AGB) stars on nuclear reaction rate uncertainties. 

For each MC iteration the reaction rate for reaction x is given 

by 

Rr(T) = cr(RS(T)~Rr,(T))+R^(T), i^cr^o (d 

Rr(T) =Cr(R^(T)-Rr(T)) + R^(T), Q>Q^-1 (2) 

where 1 ^ ^ — 1, and K^, R^, Rf are the upper, nominal, and 
lower rates respectively. We choose a Gaussian distribution to the 
random variable which peaks at zero and has an e-folding dis- 
tance at — .55 so that the extreme values are populated at 1.8cr. 
Note that while the distribution on is symmetric, the distribution 
of R(T) is not. 



2.3 Details of the computed models 

With the stellar evolution code we generated seven sequences of 
stellar structure models for a range of initial mass and metallicity. 
Each sequence was followed from the beginning of the MS to just 
after the end of the Idup. The stellar model sequences together with 
the mass coordinate of the deepest penetration of the Idup (Mdup) 
and the maximum temperature achieved at that mass coordinate are 
listed in Tablejl]. The maximum envelope temperature was achieved 
at about | of the MS lifetime for all cases. The envelopes of the 2- 
3 Mq cases receive matter from the core due to core convection 
on the early MS. For these cases. Table |l] also shows the maximum 
temperature achieved at the centre of the star (Tc) from the begin- 
ning of the MS to when the core convection no longer overlaps with 
the envelope. These temperatures are the maximum that part of the 
surface material, after Idup, was exposed to. 

Out of the full stellar evolution model sequences about 32 to 
36 models for 1-1.5 Mq cases and 82 to 1 19 models for 2-3 Mq 
were chosen for post-processing. The bulk of extra models used 
in the 2-3 Mq cases were put on the first half of the MS, when 
core convection was retreating, leaving behind material processed 
in the core. At least 4 models were put near the time of deepest 
Idup to ensure that a precise mixing depth was achieved. We post- 
processed the whole star for 2-3 Mq but only the envelope layers 
affected by the Idup for 1-1.5 Mq. In this study we ignore any 
extra mixing proces ses that might occur later on the RGB, like coo l 
bottom processing (Wasserburg, Boothroyd, & Sackmann 1995). 



ity is lower. For the lower metallicity cases we have scaled all three 
oxygen isotopes with Z, not taking in to account the anti-correlation 
of the ^''0/^*0 with Z reported by Timmes, Woosley, & Weavei 



The choice of stellar cases spans 1 Mq to 3 Mq with solar metal- 
licity and additional cases for 1 Mq and 2 Mq where the metallic- 
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(1995). To facilitate comparison, the initial oxygen isotope ratios 
were chos en to be ^^O /^'^Oj = 2465 and ^^^O/^^Oi = 442, the 
same as inElEic (1994). 

Fig. Ill shows example abundance profiles produced by the 
post-processing code using the nominal reaction rates, together 
with the Idup mark. One can see that ^^O is particularly sensitive to 
the maximum depth of the Idup for the 1 Mq and 1.5 Mq cases 



because offlieshar£_^ise_^i^ 
]r990| ). The reduced 



O at the bottom of the 



dre et al. 



TS 



ldu£_(La 



O at the bottom of the envelope. 



from where core convection retreated, enhances the sensitivity to 
the maximum depth of the Idup for the 2 Mq and 3 Mq cases. 

The seven CNO cycle reactions which are included in the anal- 
ysis are listed in TableH in decreasing order of Rn(logT = 7). 
The reaction ^''0(p,7) F is followed by a rapid f3^ decay and 
draws from a large reservoir of ^^0 to produce ^^O. Reaction 
^^0(p,7)^*F consumes '^^O to produce ^*0. Branching off from 
this series are reactions ^*0(p,a)^^N and ^*0(p, 7)^®F, which 
consume ^**0, and ^^0(p, a)^'*N, which consumes ^^O. For the 
1 Mq cases, reaction ^®0(p, 7)^'^F is slow enough and ^"^O is 
abundant enough that the ^^O reservoir is effectively not changing 
in the envelope. This is not true for 2-3 Mq (see Fig. |l|), and there- 
fore reactions ^^N(p, 7)^^0 and ^^F(p, a)^^0 should receive at- 
tention in these cases because they affect the production of ^®0. 
Fig. ^ shows the nominal, upper and lower limits on these reaction 
rates for the relevant temperature range. 



3 RESULTS 

The primary result from the MC simulations are quantified model 
uncertainties in the oxygen isotopic ratios due to reaction rate un- 
certainties. Secondly, we assess the individual contribution from 
each of the considered reaction uncertainties on isotopic ratio un- 
certainties. 

An MC simulation was carried out for each of the cases in Ta- 
ble |l] and the results are summarized in Table I and Fig. ^. The iso- 
topic ratio uncertainties are derived from the distribution of points 
in triple isotope plots. The upper error bars for ^®0/^*0 are de- 
fined to be the location of a horizontal line separates 2 per cent of 
the points above the line, and 98 per cent of the points are below. 
The lower error bar is where 98 per cent are above and 2 per cent 
are below. Error bars for ^^O/^^O are defined similarly. Therefore, 
the total number of MC points that lie outside of the error bars is 
8 per cent. This scheme helps ensure that the errors bars are 
nearly invariant to the number of iterations. There were 200 itera- 
tions for each 1-1.5 Mq case, and at least 400 iterations for each 
> 1.5 Mq case, which are already more time consuming due to 
a greater number of input models. More iterations are needed for 
the > 1.5 Mq cases because there are 3 instead of 2 influential 
reaction rate uncertainties. 



3.1 The triple isotope plots from the MC simulation 

Fig. 1^ shows a typical result from an MC simulation. Each point 
on the triple oxygen isotope plot represents the isotopic ratios af- 
ter the Idup for one iteration. The set of for each iteration al- 
lows one to locate points with reaction rates far from the nominal 
value. In Fig. ^ a symbol associated with reaction x is over-plotted 
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if > 0.7 (i.e. > lAa) is satisfied for that reaction. Multi- 
ple symbols may be over-plotted for individual iterations. In this 
way the oxygen isotopic ratios from iterations for which reaction 
X was near an upper or lower limit are flagged. The distribution of 
the flagged iterations can reveal the relative importance of the un- 
certainties of the reactions considered. Symbols for reactions that 
do not dominate the isotopic uncertainties can appear everywhere 
in the cloud of points, indicating that a rate far from the nominal 
value for these reactions does not cause any preferred location in 
the plot (reaction ^^N(p,7)^®0 for example). The lack of sym- 
bols from reaction ^*0(p, q)^^N in a horizontal band through the 
central region of the plot indicates that its reaction rate uncertainty 
overwhelmingly dominates the others in affecting the ^"O/^^O er- 
ror. Reaction ^*0(p, a)^^N causes scatter on the ordinate because 
it destroys ^*0. Fig. ^ shows this more clearly. 

For the 3 Mq case shown, reactions ^''0(p,7)^^F and 
^^0(p,ci)^*N are responsible for the scatter on the abscissa be- 
cause they most strongly affect ^^O. This fact is not clear from Fig. 
^because neither one clearly dominates. With random there are 
many iterations where there is a cancelling effect, leading to points 
in the middle part of the cloud. In order to demonstrate this we do 
15 separate runs of the post-processing code. In the first run, all 
rates are held at nominal. In the subsequent 14 runs, each of the 7 
reaction rates takes its upper, then lower value while the other rates 
remain at the nominal value. The result of the 15 runs is illustrated 
in Fig.^for the 3 Mq case. Reaction ^^0(p, a)^'*N affects the ab- 
scissa and reaction ^ 0(p, 7)^^F affects the ordinate secondly, but 
primarily affects the abscissa, as one might expect from inspection 
of the CNO network. Fig. ^ shows a similar 15 combination run 
on the 1.5 Mq case. For masses ^ 1.5 Mq the envelope does not 
receive matter processed in the core, and hence the surface abun- 
dances are a product of nuclear processing at cooler temperatures. 
As a result, uncertainty of the ^^0(p,ci)^*N reaction loses all of 
its effect on ^''O/^ ' O uncertainty for stellar mass ^1.5 Mq cases. 
Fig. shows the absence of points in the centre for the dominant 
reactions. 

The point distribution in the MC triple isotope plots of course 
depends on the choice of the e-folding distance for the random vari- 
able The distributions typically have patterns due to the asym- 
metric nature of reaction rate uncertainties. For example the bi- 
modal distrubution in Fig. n (high point density in the lower half 
and low density for larger 0/^**0) is a consequence of the large 
upper limit for reaction ^*0(p, a)^^N. 



3.2 Discussion of reaction uncertainties 

Reaction ^®0(p,7)^^F has a relatively low uncertainty (Fig.^. 
For example at log T — 7.4 the upper and lower limit according 
to the NACRE compilation are ±30 per cent. This reaction de- 
rives its effect from the large (virtually constant) reservoir of ^"O 
which leads to a strong production rate of ^^O. Although the reac- 
tion rate of ^^0(p, a)^*N is about two times greater than that of 
^^0(p,7)^^F, ^^O is produced because the abundance of ^''O is 
about 3000 times greater than that of ^^O. Hence, the uncertainty 
in reaction ^''0(p,7)^^F is more important to the surface abun- 
dance of ^^0.0 In the past, reaction rates for ^^0(p, a)^'| N and 
^^0(p, 7)^*F were very uncertain (see Sect. 3.2 |e1 Eid||l994| ). With 
the modem, smaller NACRE uncertainty, uncertainty in reaction 



^^0(p, a)"N is still relevant to the ^^O/^'^O Idup predictions for 
the 2-3 Mq cases. Modern uncertainty in reaction ^^0(p, 7)^*F, 
however, is not significant in the cases studied here. 

For the 2-3 Mq cases, reactions ^^0(p,a)^*N and 
^^0(p,7)^^F (starred triangles and triangles respectively in Fig. 
^ have a similar influence on ^^O/^^O. In the 2-3 Mq cases the 
combined effect of the uncertainty in these two reaction rates can 
cooperate to give some extreme ^"O/^^O values as well as partially 
cancel to give points in the central part of the cloud. 

Reaction ^*0(p, q)^^N strongly effects the isotopic abun- 
dances because its uncertainty is very large at the envelope temper- 
atures of these stellar cases. Also note that reaction ^*0(p, 7)^^F 
has no noticeable effect on ^*0. Even though the error for this re- 
action is very large, reaction ^**0(p, a)^^N is faster and therefore 
more important to the destruction of ^*0. 

The effect of stellar mass and metallicity on the isotope ratio 
uncertainties is detectable mostly in ^^O/^'^O. The ^®0/^®0 uncer- 
tainty is typically -|-20/ — 5 per cent due to reaction ^*0(p, a)^^N. 
For cases > 1.5 Mq, with higher temperatures (2 Mq, Z=0.001 
and 3 Mq, Z=0.02 especially) greater ^''O destruction and ^^O 
production can be seen at the surface. Hence, the destruction rate 
of ^^O rivals its production from ^^O, so the effect of uncertainty 
in reaction ^^0(p, a)^*N is visible at the surface. 

The interaction between the core and envelope convection for 
the > 1.5 Mq cases means different temperature environments 
are seen by some of the dredged up material, than are seen by the 
^ 1.5 Mq cases. This adds to the complexity and has a very no- 
ticeable impact on the propogation of reaction rate uncertainty to 
surface abundances because it means uncertainties at a wider vari- 
ety of temperatures become important. 



3.3 Comparison with observations and other models 

Our isotopic ratios, with the nominal r a tes, ag ree with results from 
Boothroyd, Sackmann, & Wasserburg ( 199^ ) that predict a strong 
dependence of ^"O/^'O on stellar mass. With increasing stellar 
mass from 1 Mq to 2 Mq, ^^0/"0 decreases from 2500 to 100 
for Z=0.02 in Fig. |^ (i.e. destruction in the envelope increases 
for > 1.5 M(i) and the ^^O production is increased in all cases). 



^0 

Our predictions show increasing ^^O/^^O with increasing mass for 
the > 2 M(7) cases, producin g a minimum at about 2 Mp), the same 
as in |E1 Eid ( 1 994 ) (see also Boothroyd, Sackmann, & Wasserburg 
1994p . Ttieinset in Fig. |8| shows model predictions by ElEidp 1994^ 



and this study (asterisks and filled boxes respectively) for 2 Mq 
and 3 Mf^^. The difference between our ^^O/^^O predictions and 



those of El Eic ( 1994 ) is partially explained by our different choice 
for the important 0(p, a)^*N reaction rate. The lower limit on 
this rate, used in this study, is similar to the rate used to predict the 
asterisks. Uncertainty in the treatment of mixing processes, which 
is not propagated with the MC simulations, likely makes the largest 
contribution to the difference. 

The dependence of oxygen abundance ratios on stellar metal- 
licity is tied to dependence of structure (temperature) to metal- 
licity. The predicted ratios for the two 1 Mq cases are indistin- 
guishable given the calculated uncertainty. The predictions for three 
2 Mq cases are distingushable. Lower metallicity produces higher 
temperatures and greater overall ^*0 destruction, causing higher 

Spectroscopic obse r vation s of the isotope ratios in RGB stars 



^ Production (or destruction) rate is r = NaNi,{av), where Na and A^^ 
are the molar densities of the reactants and (uv) = Ris the reaction rate. 



from [Harris & Lamberj ( |1984[ ); [H arris, Lambert, & Smith] ( |1988| ) 
are also shown in Fig. Q and Table ^. These stars are either ascend- 
ing the RGB for the first time or are undergoing core He-burning 
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on the blue loops. Hence, the surface abundances of these stars 
have been changed by the Idup. The mass-^^' O/^^O relatio nship 
is roughly obeyed by these observatations (see El Eid 1994, for a 
discussion). 

The prediction band in Fig. ^ was from an initial ^"O/^^O 
which is 10 per cent smaller than the solar value, giving pre- 
dictions that overlap with all but one of the spectroscopic obser- 
vations. However, the number of post-ldup observations below 
solar (^''0/^*0=498) and below the band predicted from initial 
16q^18q_442 is worth noticing. When there is a slight shift in the 
initial oxygen isotopic ratios, there is a proportional shift in the 
modeled p ost-ldup ratios (as shown by Boothroyd, Sackmann, & 
Wasserburg |1994| ). Hence, if the model prediction of about a 15 
per cent increase in ^''0/^**0 after Idup is correct, then some of 
the stars with spectroscopic data in Fig. |^ (circles) must have had 
initial ^''0/^**0 values ~ 30 per cent less than solar (super-solar 
metallicity). In addition, the post-ldup spectroscopic data shown 
in Fig. P roug hly correlate with [Fe/H] measurments compiled by 
Tayloi (1999) (Table but indicate a pre-solar metallicity. There 
is currently no reason to believe that is not depleted in the en- 
velope during the MS. A possible explaination is additional mixing 
above the core that further depletes ^^O in the envelope, leading to 
an decreased ^^0/^*0 after Idup. Testing this hypothesis will be 
appropriate when more precise spectroscopic observations become 
available. 

Pre-solar meteoritic inclusions have been linked t o giant stars 



due t o their specifi c isotopic abundance signatures (Huss et. al 



1994; Zinner 



19981. The rather large spread of their oxygen iso- 



topic ratios implies additional mechanisms at work (like intital iso- 
tope ratio variations or extra mixing processes). The variations of 
isotope ratios in grains are much larger than the model uncertainties 
due to nuclear reaction rates. However, better nuclear data would be 
required to improve the identification of these extra mechansims. 



4 CONCLUSIONS 

In this paper we quantify the uncertainty in oxygen isotopic ratios 
due to uncertain reaction rates. These results may help to moti- 
vate and prioritize new laboratory measuremtents of reaction rates. 
Reaction rate uncertainty for ^''0(p, 7)^^F and ^**0(p, q)^^N are 
significant to oxygen isotope ratios for all of the stellar cases stud- 
ies here. Reaction ^^0(p, q)^''N competes with ^''0(p,7)^^F for 
dominance in ^^OI^^O uncertainty for the 2-3 Mq cases because 
of the interaction between core mixing and subsequent envelope 
mixing. Of course, this result is based on the assumption that the es- 
timated uncertaities in the reaction rates (NACRE upper and lower 
limits) are appropriate. 

In general, reactions with large uncertainties, and reactions 
with slow reaction rates but large production/destruction rates and 
any uncertainty create the largest uncertainty in isotopic ratios. The 
MC scheme demonstrated here efficiently finds the most problem- 
atic rates and provides a means of quantifying uncertainties for a 
particular stellar environment. 

The model uncertainties calculated here are from reaction rate 
uncertainties only. The total model uncertainty in oxygen istopic 
ratios definitly has contributions from uncertainty in mixing pro- 
cesses. This problem will become more tractable as the uncertain- 
ties due to reaction rates are made smaller. 
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Table 1. The seven different stellar cases considered with the maximum temperature achieved in the envelope during the MS. The mass coordinate of the 
deepest penetration of the Idup is M^up and Tc is the central temperature. 



mass 


metallicity 


Mdup 


max[logT(Mdup)] 


max[logTc] 


(Mq) 




(Mq) 


[K) 


[K) 


1.0 


0.01 


0.23 


i.n 




1.0 


0.02 


0.225 


7.11 




1.5 


0.02 


0.25 


7.22 




2.0 


0.001 


0.39 


7.36 


IM 


2.0 


0.01 


0.32 


7.32 


7.35 


2.0 


0.02 


0.31 


7.29 


7.29 


3.0 


0.02 


0.49 


7.36 


7.39 



Table 2. Reactions included in the MC scheme, with nominal, lower and upper limits for reaction rates at log T = 7.0 and log T = 7.4. 



log{iJ„(7.0)) log (R„ (7.4)) 
{sec~^mol cm~^) 



i**0(p, o)15n 


-19.98: 


15n(p,7)1'50 


-20.36: 


180(p,7)19F 


-23.00: 


"F(p,a)i'5o 


-23.75: 


"0(p,a)i'*N 


-23.77: 


"0(p,7)"F 


-24.17: 


i60(p,7)"F 


-24.16: 



+ 1.13 _ -I -I 0-1+0.30 

-0.70 ^^•"^-0.26 

+0.07 _ -.r, q-,+0.09 

-0.09 ^^-^^-0.12 

+2-42 „„+0.17 

-0.78 -^^•"^-0.04 

-0.16 -^^-^^-0.16 

-0.15 ^'-'•""-0.09 

+0.08 _ 1S+0.12 

-0.10 ^'^•^°-0.12 
-l^/tl+O-ll 

-0.15 ^^-^^-0.15 



Table 3. The modeled surface oxygen isotopic ratios after Idup, with errors from the MC simulations. 



mass metallicity, Z ^^O/^^O ^'^O/^'^O 
(Mq) 



1.0 


0.01 


2405^17 


470+t^ 


1.0 


0.02 


2410116 


469+t^ 


1.5 


0.02 


12601^3° 




2.0 


0.001 






2.0 


0.01 




575ir 


2.0 


0.02 


1 1 5+3** 

ilO_2o 


565+^^ 


3.0 


0.02 


268l|» 




initial 




2465 


442 
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4. Selected observations of oxygen isotope ratios. The first 8 are measured from spectroscopic observations of RGB stars. The two 
meteoric observations were extracted from individual dust grains that formed around AGB stars. 



Object massXMQ) "*0/^^0 ^''O/^^O ref'^ [Fe/Hf 



a Ah 


1 


J^O_j20 


450 +'*° 
^-'^ -no 


1 


-0.221 


13 Gem 


1.3 


240 til 


•"^^ -no 


1 


-0.003 


a Ser 


1.4 




400 


1 


+0.098 


/3UMi 


1.5 


510 til 


4401^° 


1 


-0.132 


a Tau 


1.5 


eootlTo 


47^ +125 
-200 


2 


-0.102 


13 Peg 


1.7 


1050 till 


475 +75 

-150 


2 




7 Dra 


2 


300 tllo 


47^ +125 
^'J -200 


2 


-0.178 


fj, Gem 


2 


+75 
J^J -150 


47^ +125 
^'J -200 


2 




/3 And 


2.5 


iJJ -50 


425 tllo 
600 11^^ 


2 




a UMa 


3 


330 trl 


1 


-0.128 


SEAL203 




355 


498 


3 




SEAL235 




472 


714 


3 




SEAL261 




337 


1190 


3 




solar 


1 


2622 


498 


4 






■ Masses have an undetermined uncertainty. 



^ Observations 1) at 5/im by Harris, Lambert, & Smith (1988) and 
2) Harris & LamberjJl984^^_^2j™teori^^ Huss, & 



Wasserburg ( 1998| ), 4) solar ratios by [Anders & Grevesse| ( 1989 1 



From Tayloi (1999) 



© 2002 RAS, MNRAS 000, 



Oxygen isotopic ratios in first dredge-up red giant stars and nuclear reaction rate uncertainties revisited 





\ 1 
V ! 


© 




16(-) 

: no / 


\ 1 

\ 1 ■ 




- 












j\ 

\ / ' 
\ / ' 

V i 
1 







\ 

\ - 
\ - 
\ - 



0.2 0.4 0.6 0.8 






1 






■ — 


















, 1 




, " r 


0.5 


1 

M/M^ 


1.5 


2 








' 1 . 










■ '-'0--,-''' 

1 / 














1 

/ 

\ \ 






"^X : 
' 1 " 


1 




2 


3 



- 7.5 
7 



Figure 1. Mass fraction profiles of H (dotted), ^^0, ^''O, and ^^O (solid), 
as well as temperature (long dash) for the 1 Mq (top left), 1.5 Mq (top 
right), 2 Mq (bottom left), and 3 Mq (bottom right) cases with Z=0.02 
near the end of the core H-burning. Also shown is the deepest penetration 
of the first dredge-up (dash-dot). 
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-14 



-16 



-20 



-22 - 
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Figure 2. The nominal, upper and lower limits of the rates for the seven 
reactions. The order of reations in the legend is the same as Table |^. 
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Figure 4. Oxygen isotope ratios for the 3 Mq Z=0.02 case when one re- 
action took the upper, lower, and nominal (circled) reaction rate while the 
other rates remained at nominal. Results from ^^F(p, a)^^0 (not shown) 
are the smallest. 
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Figure 5. Same as Fig. ^ for the 1.5 Mq Z=0.02 case. 



Figure 3. The surface oxygen isotope ratios for 3 Mq and Z=0.02 from 
an MC simulation with all seven reactions and 500 iterations. All points 
are plotted in small filled squares and over-plotted are points from iter- 
ations where \^f\ > 0.7 for each reaction; ^^N(p,7)^^0 (diamonds), 
^^0(p,7)^'^F (open triangles), ^'^0(p, q:)^'*N (starred triangles up), 
l^O{p,7)l**F (asterisk), ^^OCp, q:)1'5N (open squares), l80(p,7)l''F 
(starred triangles down), ^^F(p, a)^^0 (starred diamonds). Note the lack 
of open tiiangles and squares in the centre (error dominant reactions 
^^0{p, 7)i^F and i**0(p, a)i^N). 



This paper has been typeset from a Tj^X/ W^Si file prepared by the 
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Figure 6. Same as Fig. g for the 3 Mq case, except only symbols from it- 

erations where | > 0.7 (open triangles), m-' / | ^ 

0.7 (starred triangles up), and '^'P'"' ^| > 0.7 (open squares) are 
over plotted. A filled circle marks the point where the nominal rates were 
used for all reactions. 
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Figure 8. Modeled isotopic ratios with propagated errors derived from MC 
simulations for the seven stellar cases (filled squares, labeled with their 
mass and metallicity). The x marks the initi al oxygen ratios be f ore Id up for 
all cases. Spectrosco pic data (circles) from [Harris & l.amherl i 1 984[): H ar- 
lis, Lambert, & Smith ( 198j ). Meteoric data from SEAL203 (open square), 
SEAL235 (open triangle) (Choi, Huss, & Wasserburo 1998 1, and solar (di- 



amond). Dotted lines mark a band where Z=0.02 stars are predicted to exist 
with our choice of initial abundances. Inset show s our 2-3 model pre- 
dictions in more detail, as well as predictions of El Eid| (1994 ) for Z=0.02, 
2 Mq (left asterisk) and 3 Mq (right asterisk). 
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Figure 7. Same as Fig. ^ for the 1.5 Mq case with 200 iterations and 
only symbols from iterations where 0(p,7) ^| ;> g. 7 (open triangles) 
and 0{P'") N| ^ g y (-Qpgjj squares) are over plotted. A filled circle 
marks the point where the nominal rates were used for all reactions. 
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